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We present a limit on the brancliing fraction for the decay S" — > 77 using data collected at the 
T(45) resonance with the BABAR detector at the PEP-II asymmetric-energy e'^e~ collider. Based on 
the observation of one event in the signal region, out of a sample of 21.3 x 10 e^e~ — > T(4S') -^ BB 
decays, we establish an upper limit on the branching fraction of 3{B^ -^ 77) < 1.7 x 10"'' at the 
90% confidence level. This result substantially improves upon existing limits. 

PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 



In the Standard Model the decay B^ -^ 77 proceeds 
via a second-order weak transition, including gluonic 
penguins, followed by annihilation (Fig. ^. Standard 
Model predictions for the branching fraction of these ef- 
fective flavor-changing weak neutral current processes 
range from 0.1 to 2.3 x IQ-^ |J|. 

Physics beyond the Standard Model can enhance this 
branching fraction by as much as two orders of magm- 
tude, particularly in the case of two-Higgs models [2|- 
Other particles from the supersymmetric spectrum can 
further modify the Standard Model expectation []| . The 
current best limit on the branching fraction for i?" ^77, 
from the L3 experiment at the CERN LEP collider, 
is B{B° -^ 77) < 3.9 X 10-'^ (90% confidence level). 

In this Letter we present an analysis based on data 
taken with the BABAR detector \^ , which operates at the 
PEP-II asymmetric-energy e+e" collider at the Stanford 
Linear Accelerator Center p|. The sample consists of 
19.4 fb~ taken at the T{4S) resonance, corresponding to 
21.3 X 10^ e+e" -^ r(4S') -^ BB events. An additional 
sample of 2.2 fb"^ accumulated 40 MeV below the T{4:S) 
resonance is used to estimate non-BB background. 

Charge conjugation invariance is assumed for all chan- 
nels quoted in this paper, and the charge conjugate reac- 
tions are included in the analysis. Quantities evaluated 
in the T{AS) rest frame are denoted by an asterisk; 6.17., 
E^ is the energy of the e"*" and e" beams in the T(4S') 
rest frame. 

The BABAR detector, a general purpose solenoidal 
magnetic spectrometer, is described in detail elsewhere 
y . A silicon vertex detector and a cylindrical drift cham- 
ber in a 1.5-T solenoidal magnetic field are used to mea- 
sure momenta and ionization energy loss of charged par- 
ticles. Electrons and photons are identified by a Csl elec- 
tromagnetic calorimeter (EMC). 

This analysis exploits in particular the information 
provided by the EMC consisting of 6580 Csl crystals, 
covering 90% of 47r in the T(45) rest frame. The en- 
ergy resolution has been measured directly with a ra- 
dioactive source at low energy and with electrons from 
Bhabha scattering at high energy. The mass resolution 
of tt'^ and rj candidates in which the two photons in the 
decay have approximately equal energy can be used to 
infer the energy resolution at an energy less than 1 GeV; 
the decay Xci -^ J I'^l provides an additional measure- 
ment at 500 MeV. A fit to the energy dependence results 
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FIG. 1: Examples of possible diagrams responsible for the 
decay B° — > 77. In these diagrams q = u, c, or t, and H is a 
hypothetical charged non-Standard Model Higgs boson. 



in geIE = (2.3 ± Q.i)%/^^fE © (1.9 ± 0.1)% 0. 

Energy deposits in the EMC are reconstructed by 
grouping adjacent crystals with energy deposits greater 
than 1 MeV into clusters. Clusters with more than one 
local energy maximum are then split into bumps. The 
energy of each crystal is divided among the bumps by 
an iterative adjustment of the centers and energies of 
the bumps assuming electromagnetic shower shapes |^. 
Next, all tracks reconstructed in the tracking volume are 
extrapolated to the EMC entrance and a track-bump 
matching probability is calculated for each pair. 

All bumps with a matching probability smaller than 
10~^ are treated as photon candidates. Photons are se- 
lected by requiring the bump shape to be compatible with 
an electromagnetic shower, and by requiring the bump 
to have a minimum energy of 30 MeV. In addition we 
accept only photon candidates which are isolated from 
any other bump in the event. This requirement selects 
against background from high-energy tt" mesons, where 
the two photons from the decay of the tt*^ meson strike 
the calorimeter in close proximity [merged tt*^). 

The BABAR detector is simulated by a GEANT-based 
Monte Carlo procedure |g that includes beam-related 
background by mixing random trigger events into the 
Monte Carlo generated events. The simulated events are 
processed in the same manner as data. The simulation 
is used to study background and optimize selection cri- 
teria, but only enters the analysis directly through the 
calculation of the signal efficiency. 

In order to select BB events, we require at least three 
tracks of good quality in the event. The quality require- 
ments for these tracks include a small impact parameter 
with respect to the collision point along the beam direc- 
tion (10cm) and transverse to it (1.5cm), a minimum 



number of 13 hits in the drift chamber and a momentum 
of p < 10 GeV/c in the laboratory frame. To help re- 
ject continuum background, the ratio of the second Fox- 
Wolfram moment to the zeroth Fox- Wolfram moment [3 
must be less than 0.9. We further require that there be 
two high-energy photon candidates with an energy in the 
T{4S) rest frame between 1.5 and 3.5 GeV. At this point, 
all remaining pairs of photons are considered candidates 
for the decay B° -^ 77. If the event contains more than 
one such B candidate all of them are kept for further 
analysis. 

After this preselection additional requirements are im- 
posed on the i?° -^ 77 candidates. Photon bumps from 
the B candidate must not contain noisy crystals or crys- 
tals which produce no signals. The second moment of the 
energy distribution around the cluster's centroid must be 
smaller than 0.002. This value has been optimized to re- 
ject remaining background from merged tt*' mesons. 

Since B mesons at the T(4S') resonance are produced 
nearly at rest, the decay B'^ — > 77 will contain two 
nearly back-to-back photons with E* « 2.6 GeV in the 
T{4:S) rest frame. This represents a clean signature 
and makes this channel relatively easy to study exper- 
imentally. We exploit this feature by considering only 
50 _^ ^^ candidates which have at least one photon with 
2.3<£;; < 3.0 GeV. 

In order to reject photons from Tr'^irj) decays we com- 
bine each photon from the B candidate with all the 
other photons in the event having energy greater than 
50(250) MeV. The resulting Tr'^{ri) candidates are re- 
quired to have an invariant mass beyond three standard 
deviations, or 3 x 8.8(18) MeV/c^, of the nominal 7r°(77) 
mass (81 . 

Reconstruction of exclusive final states from B mesons 
produced at the T{AS) resonance benefits from the beam 
energy constraint E^ = E^. Thus, in the T(4S') rest 
frame the energies of the B meson decay products must 
add up to the beam energy. We calculate the energy 
difference AE = E, 



7,1 



^;,2 



E^ between the candi- 



date -8° meson and the beam energy in the T(4S') rest 
frame. The distribution of this quantity peaks at GeV 
for true B mesons, and has a tail towards negative AE 
due to shower leakage in the EMC. The resolution in AE 
is obtained from signal Monte Carlo events with a fit of 
the AE distribution to an empirical function [Qj and is 
(TAB = 73 MeV. 

The B meson mass resolution is improved with the use 
of the beam energy constraint. We use the beam energy 

substituted mass toes = wi?^^ — (p* ^ -I- p* 2)^- The 
resolution on ttt-es is obtained from signal Monte Carlo 
events with a fit of the tties distribution to an empirical 
function p and is CTtoes ~ 3.9 MeV/ c^. 

For the purpose of determining numbers of events and 
efficiencies a rectangular signal region is defined. This 
region extends 2a in AE about OMeV and 2a in tties 



TABLE L Cumulative signal reconstruction efficiency as se- 
lection criteria are applied. The first row shows the cumula- 
tive event selection efficiency. The additional rows give indi- 
vidual contributions to the B candidate selection efficiency. 
The cumulative signal reconstruction efficiency is given by 
the product of event selection and final B reconstruction effi- 
ciency. 



selection criteria 



efficiency [%] 



cumulative event preselection 



39.8 



photon energy E* 92.9 

bump quality and second moment 86.8 

■K° and ri veto 72.8 

|cos(6Ib)| and |cos(6It)| 40.1 

signal region 27.0 



cumulative 



10.7 



about the nominal mass tti^o of the B'^ meson. 

The search for S" — + 77 was performed as a blind 
analysis by hiding a 3a region in AE and toes in on- 
resonance data until the development of the selection 
procedure was complete. This allows optimization of the 
selection and estimation of the background without the 
bias of knowing the number of events in the signal region. 

Monte Carlo studies indicate that the main back- 



ground arises from the process e^e —^qq {q 



d, s), 



referred to as continuum background and modeled with 
the JETSET event generator |[lOj. Such events exhibit a 
two-jet structure and contain high momentum, approx- 
imately back-to-back tracks. One source of background 
includes photons from initial-state radiation, others are 
photons from 7r° —^ 77 and rj — > 77 decays, where the de- 
cay is very asymmetric in the final-state photon energy. 
Background from merged tt'^ mesons is negligible. 

To reduce continuum background, we calculate the an- 
gle 9^ between one of the photons (chosen randomly) 
of the B'^ candidate and the thrust axis of the remain- 
ing tracks and neutral bumps in the event. The dis- 
tribution of I cos 6^ I is uniform for signal events and 
strongly peaked at 1 for continuum background events. 
We also calculate the angle 9'^ between the momentum 
vector of the i?" candidate and the beam axis in the 
T{4:S) rest frame. The distribution of |cos0^| is uni- 
form for continuum background and follows a sin 9'^ 
distribution for signal events. The requirements for both 
I cos 9^ I and | cos 9g \ have been optimized to maximize 
the statistical significance Ns/y^Ns + Nb, where Ns is 
the number of signal candidates expected assuming for 
the branching fraction B{B^ -^ 77) = 1 x 10~^ [l| 
and Nb is the expected number of background candi- 
dates determined from continuum Monte Carlo simula- 
tion and off-resonance data. We require | cos9^\ < 0.57 
and |cos6'|j| < 0.81. If more than one B meson candi- 
date per event remains after this selection, which occurs 
in less than 0.1% of the events analyzed, we select the 
candidate with the smallest \AE\. After all these selec- 
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FIG. 2: Energy diflFerence AiJ between the candidate B 
meson and the beam energy in the T(45') rest frame versus 
beam energy substituted mass rriES for on-resonance data. 
We observe one event in the signal region, outlined as a black 
dashed box about AE — GeV, consistent with the expected 
background. The dashed box on the left shows the sideband 
used for background estimation. 



tion criteria the overall efficiency for B'^ -^ 77 decays 
is determined from the Monte Carlo simulation to be 
(10.7 ±0.2)%, where the error is purely statistical. Table 
^shovifs the cumulative signal reconstruction efficiency as 
the selection criteria are applied. 

A single event in the on-resonance data meets these 
selection criteria, as shovsfn in Fig. |2| A number of 
exclusive decay modes that can mimic B^ ~^ 77 de- 
cays have been studied with high statistics (equivalent 
to 1.2 — 1.7 X 10"'fb~ assuming branching fractions of 
the order of 10"®). We expect negligible contributions 
from B^ —> Tjrj, K*^^, P^li and tt'^tt'^, and a combined 
contribution of 0.7 x 10~^ events from B^ -^ p^{t:^'it^)^ 
and B^ -^ u;(7r''7)7. To further explore the question of 
remaining background in the signal region, we define the 
grand sideband consisting of a rectangular region within 
the limits -1.0 < A^ < 1.0 GeV and 5.20 < rriEs < 
5.26 GeV/ c^ (see Fig. |21 left dashed box). In this region 
we find a prediction of 34 ± 9 events from continuum 
Monte Carlo simulations, in good agreement with the 
observation of 43 ± 7 (44± 20) events from on-resonance 
data (off-resonance data of 2.2 fb~^ scaled to the full 
analyzed luminosity of 19.4 fb~^). We parameterize the 
background using on-resonance data. The background in 
AE is parameterized in the grand sideband with a first- 
order polynomial (see Fig. |3^) ; the background in rriEs is 
parameterized in the lower sideband, which is a rectangu- 
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FIG. 3: a) Fit to the AE distribution in the grand sideband 
to a first-order polynomial; b) fit of the rriES distribution in 
the lower sideband with the ARGUS function (iJl . See text 
for the definition of the sidebands. 



lar region within the limits — 1.0 < AE < — 0.2 GeV and 
5.20 < rriES < 5.29 GeV/ c^, with an empirical threshold 
function first employed by the ARGUS collaboration [llj 
(see Fig. Et). Both parameterizations describe the cor- 
responding distribution very well with a x^, normalized 
to the number of degrees of freedom, of about 0.8. Us- 
ing this parameterization we are able to extrapolate the 
on-resonance grand sideband data into the signal region 
and find an expectation of 0.9lg'3 events. This is consis- 
tent with the hypothesis that the observed event in the 
signal region is due to continuum background. Never- 
theless, we choose to quote a conservative upper limit, 
assuming that the observed event in the signal region is 
in fact due to the decay B^ -^ 77. We use Poisson statis- 
tics to set an upper limit on the branching fraction. The 
upper limit on the branching fraction B is obtained from 
B = NuL/{e ■ [Ngo + Ngo)), where Nul is the upper 
limit on the number of observed events, e the signal re- 
construction efficiency of (10.7 ± 0.2)% and Ngo+Nao is 



BO 



the number of produced i?° and B^ mesons. Ngo + N 
is equal to the number of T(45) events since we assume 
the number of i?° i3° events to be 50 % of the number 
of produced T(4S') events. This yields an upper limit 



TABLE II: Summary of systematic uncertainties on the signal 
efficiency and the number of produced T(4S') as an error on 
the branching fraction determination. The total systematic 
uncertainty is the sum of the individual contributions added 
in quadrature. 



systematic uncertainty 



(Ag/g)% 



Number of produced T(4S') 1.6 

Photon detection efficiency 6.5 

rj veto 2.0 

-!T° veto 2.0 

AE selection 5.3 

rriEs selection 2.6 

Track finding efficiency 1.8 

Number of signal Monte Carlo events 2.0 



Total 



9.6 



at the 90% confidence level. This result improves the 
existing limit [j] by over a factor of 20. 
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on the branching fraction, based on statistics alone, of 
B{B" ^ 77) < 1.7 X 10-6 at the 90% confidence level. 

Systematic effects arise from the modeling of the signal 
efficiency and the estimation of the number of B mesons 
in the data sample. A summary of all systematic errors 
is provided in Table ITTl The most significant sources are 
the photon detection efficiency and the AE selection due 
to the uncertainty in the photon energy scale and pho- 
ton energy resolution. The systematic uncertainty on the 
photon detection efficiency has been determined from a 
study which compares the precisely know ratio jS] of the 
T -^ nn^Vr and r -^ -K-K^n^Vr rate in Monte Carlo events 
and data. This uncertainty depends on the event mul- 
tiplicity, whose effect is estimated by embedding photon 
bumps from radiative Bhabha events into both generic 
B meson and generic B meson Monte Carlo events. The 
uncertainty in the energy scale is estimated with a study 
of symmetric 77 ^ 77 decays, where both photons are 
within a narrow energy range. Systematic shifts of the 
reconstructed t] mass from the nominal value measure the 
uncertainty in the energy scale in this energy range. 

In order to include our systematic uncertainty in the 
determination of the upper limit, we follow a prescription 
given by [I^l- The branching fraction B is calculated as 
B = n/S, where n is the number of observed events and 
S — 2.3 X 10^ is the sensitivity, given by the product 
of the number of B^B^ events and the overall B^ -^ 
77 selection efficiency. Assuming a normal distribution 
for the uncertainty in 1/S, the systematic uncertainty 
is accounted for by convoluting the Poisson probability 
distribution for the assumed branching fraction with a 
Gaussian error distribution for 1/ S. Our total systematic 
uncertainty of 9.6% included in this way has a negligible 
effect on the upper Hmit. 

In summary, we have performed a search for the decay 
50 -^ ryj_ •y^e observe one event in the signal region and 
infer an upper limit on the branching fraction of 

S(B°^77) < 1.7 X 10^6 
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